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Abstract. We study graphene with an adsorbed spin texture, where the localized spins 
create a periodic magnetic flux. The latter produces gaps in the graphene spectrum and 

O , breaks the valley symmetry. The resulting effective electronic model, which is similar 

to Haldane's periodic flux model, allows us to tune the gap of one valley independently 
from that of the other valley. This leads to the formation of two Hall plateaux and a 
quantum Hall transition. We discuss the density of states, optical longitudinal and Hall 
conductivities for nonzero frequencies and nonzero temperatures. A robust logarithmic 

"j^ . singularity appears in the Hall conductivity when the frequency of the external field 

agrees with the width of the gap. 
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1. Introduction 

Transport properties of neutral graphene are characterized by a semimetallic behavior 
with a point-like Fermi surface at two valleys. This can be changed by external electric 
fields in gated graphene to more classical transport with a circular Fermi surface or by an 
external magnetic field that opens gaps in the spectrum [H [2] . The latter leads to the 
formation of Hall plateaux in the Hall conductivity. Another possibility is to introduce a 
gap either by chemical doping (e.g. with hydrogen [3]) or in bilayer graphene with a double 
gate |11 [U El [3 El [12]. Chemical doping is an interesting direction for modifying graphene 
because it leads to a rich field with new properties. Although a standard technique for 
semiconductors, it has been applied to graphene only recently P |9l [101 [TTJ, [121 [131 HU [15]. 
As a special case of chemical doping one can use atoms that carry a spin to create a spin 
texture on graphene. These localized spins alter the transport properties of graphene 
significantly. This will be discussed in this paper. 

The quantum Hall effect in graphene was observed in a number of experiments, 
cf. Ref [H [2]. It is usually associated with the presence of an external homogeneous 
magnetic field which separates the electronic spectrum into Landau levels. Without an 
external magnetic field each valley of the honeycomb spectrum provides a quantum Hall 
step [17]. However, the system is time-reversal invariant. This implies that the Hall 
conductivity vanishes, since the contribution of the two valleys to the Hall conductivity 
cancel each other. A magnetic flux, however, can cure this problem by breaking the time- 
reversal invariance. A possible way to observe the quantum Hall effect is by introducing a 
periodic magnetic flux. This was suggested by Haldane [16]: A staggered magnetic flux, 
in combination with nearest and next-nearest neighbor hopping and a staggered potential 
on the honeycomb lattice, creates two Hall plateaux with Hall conductivities a xy = ±e 2 /h. 
The periodic magnetic field leads to a staggered flux with zero net flux in each unit cell. 
It affects the hopping matrix elements by creating a phase factor which enables us to 
change the signs of the Hall conductivities of each node independently by changing the 
magnetic field strength. The effects of the staggered potential and the staggered magnetic 
flux can also be understood in terms of symmetry breaking: The Brillouin zone has a six- 
fold energetic degeneracy due to the vanishing energy at the corners of its hexagonal 
structure. A staggered potential breaks the inversion symmetry [171 [16] but preserves the 
six-fold energetic degeneracy, because all six corners acquire the same gap. The staggered 
magnetic flux, however, reduces the six-fold degeneracy to a three-fold degeneracy (i.e. 
it provides inter-valley symmetry breaking) because it affects the gap of the two Dirac 
nodes differently. 

The importance of sublattice and inter-valley symmetry breaking in graphene for 
applications has been widely recognized. Recently the effect of the broken inversion 
symmetry in graphene was addressed to in context of topological [18] and anomalous 
thermoelectric [19] transport, as well as valley-dependent optoelectronics [20] . 

In this paper we propose a model that is based on a honeycomb lattice, where one 
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Figure 1. Spinor atoms on a unit cell of graphene. Tilted spin textures can provide an 
effective staggered magnetic flux [24j [25j [26] . 



sublattice is occupied by localized spins. This leads to an effective hopping parameter 
with a Berry phase and reflects the equivalence of the spin texture with an effective 
periodic magnetic flux in the tight-binding model. It is similar but slightly different 
from Haldane's periodic flux model. The inter-valley symmetry is lifted by introducing 
only a next-nearest neighbor hopping on one sublattice and a staggered magnetic flux. 
This choice is motivated by the fact that graphene on a substrate can adsorb atoms 
only on one side. In this case the adsorbed atoms occupy only next-nearest neighbor 
sites on the honeycomb lattice, i.e. they occupy one sublattice. Such a possibility was 
experimentally observed for hydrogen on graphene [3]. As a result, the adsorbed atoms 
modify the overlap integrals on their sublattice and may even create an additional next- 
nearest neighbor hopping. Moreover, if the adsorbed atoms have a magnetic moment, 
these moments can form a localized spin texture due to double exchange [211 1221 [23] . It is 
known that such spin structures provide an effective Berry phase for the electron hopping 
rate J2U [251 I2S]. The latter has the same effect as a periodic (staggered) flux, where the 
flux depends on the tilting angle of the localized spins (Figs. HJ [2]). 

2. Model: electronic hopping on a spin texture 

We consider localized S = 1/2 spins on a graphene lattice, where the latter is assumed to 
be flat. As a reference system we choose a basis of eigenspins oriented in ^-direction with 
the lattice in the x-y plane. In other words, the z component of the localized spin 1/2 
operator satisfies S z \±) = ±|±). In terms of Pauli matrices a = (a x , a y , a z ) this reads 

a z s± = ±s± (1) 

with the two-component spinor basis {s + ,s_}. In general, a spin state can have 
a local orientation, such that it is an eigenstate to fij ■ a with the 3D unit vector 
fij = (sin 8 cos (j), sin 8 sin (j), cos 8). The matrix of the spin operator fij ■ a reads in the 
z-oriented basis 

/ cos 8j e l ^ j sin 9j \ 



J \ e 1 ^ sin 9i — cos 8 



j 



and 



Sj = nj ■ <js + = e" J 



cosdj 



sin 8j 
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Here the tilting angle 8j and the x-y rotation angle <j)j refer to the change of the 
quantization direction, relative to the spin state s+. 

Then the electronic Hamiltonian on a honeycomb lattice must describe the interaction 
of the electronic spin with the localized spins. This leads to a spin-dependent hopping 
amplitude |2H 123 I2S] : 

H = J2 H n c \ c j + h.c. (2) 

with 

H tj = tBi ■ Sj = te ia ^ cos(%/2) (3) 

and electronic creation (annihilation) operators c\ (cj). The Berry phase ay is given 
as [261: 



5^(^/2)5^(^/2)5^(0,-0,) 

sma v~ cos(%/2) [V 

which implies ay = —dji. (It should be noticed that the relation (j4j) means that ay is the 
solid angle spanned by s,, Sj and the z axis [26] .) The phase ay vanishes for vanishing 
tilting angles 8i, 8j (i.e. for ferro- or antiferromagnetic states) and for 0, = 0j (i.e. when 
the spin projections on the x-y plane are parallel. Eq. (|3]) gives us the general expression 
for the hopping of the electrons in graphene with an additional spin texture, where the 
latter is characterized by the tilting angles and x-y rotation angles. If we assume a 
specific configuration of localized spins, the corresponding quantization axis of the spin Sj 
is fixed by the angles 8j and <pj. A special case, where all angles Oij between next-nearest 
neighbor spins and all tilting angles 6j are equal (except for an irrelevant sign change of 
Oij for different pairs i,j) is depicted in Fig. [2j This is similar to the situation on the 
Kagome lattice in Ref. [26]. Moreover, the angles 0j are rotating by 2ir/3 on the lattice 
such that their differences in the expression of the phase ay in Eq. ()1]) are ±2ir/3 for next 
nearest neighbors. This implies a global renormalization of the hopping rate t — > t cos 8^ 
and a Berry phase that changes only its sign from lattice bond to lattice bond. 

2.1. Tight-binding model on the honeycomb lattice 

We begin with the tight-binding Hamiltonian of monolayer graphene (i.e. for a honeycomb 
lattice) and return to the effect of the spin texture later. Then in Fourier space we can 
write for the Hamiltonian matrix in sublattice representation 

H = h 1 a 1 + h 2 a 2 + h 3 a 3 , (5) 

where the off-diagonal Pauli matrices describe the hopping between the two sublattices. 
The diagonal Pauli matrix term with h 3 describes processes on the same sublattice that 
can include next-nearest neighbor tunneling. The eigenvalues of the Hamiltonian are 



E ± = ±E k = ±Jhj + hl + hi (6) 
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Figure 2. Spin texture on the triangular sublattice with fixed tilting angles (all x-y 
projected spins have the same length) and x-y rotation angels with multiples of 2n/3. 
The + (-) signs indicate a positive (negative) magnetic flux through the corresponding 
(up or down) triangle, whose value is determined by Eq. (j4|. 

Specifically for the honeycomb lattice we have for the nearest-neighbor terms 

3 3 

K = -£^cos(a r k), /i 2 = -t ^ sin(a- • k) (7) 

i=l i=l 

with the basis vectors of the honeycomb lattice a^ given by 

a 1 = a(0,-l), a 2:3 = | (±V3, l) , (8) 

where a denotes the lattice constant. In absence of a gap, i.e. for h 3 = 0, the spectrum 
defined in Eq. ([6]) vanishes at nodal points whose positions in the Fourier space are given 
by the vectors 

bf = — ?U(±1,0), b5 = —^(-l,±>/3), bf = —?=-(l,±V3).(9) 

Assuming a uniform gap, i.e. h 3 = m, the Hamiltonian defined in Eq. (jSJ) can be 
approximated in the vicinity of the nodal points by the Hamiltonian of a massive Dirac 
fermion 

H ps vk 1 cr 1 + vk 2 a 2 + rncr 3 , (10) 

where v = \/3ta/2 = 1 denotes the Fermi velocity of electrons. With this the eigenvalues 
read 



E ± = ±E k PS ±Vk 2 + m 2 . (11) 

Now we include the spin texture, described by the modified hopping term of Eq. (|3j), 
by the phase variable r . r / which is the Berry phase of Eq. (00). Using a spin texture 
as shown in Fig. [21 we get a Berry phase <p which is the same for all three next-nearest 
neighbor vectors Cj (j = 1,2,3 (cf. Fig. [3]) 

Cl = V3a (1, 0) , c 2j3 = ^ (-1, ±V3) , (12) 
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Figure 3. Magnetic flux through the graphene unit cell corresponding to the potential 
Eq. (|14[) . Arrows indicate second- neighbor hopping processes and <po ~ h/e is the flux 
quantum. 



because the tilting angle is fixed and the x-y rotation angle of the spins is ±27r/3 between 
next-nearest neighbor sites. Then the next-nearest neighbor hopping term reads 



Xr 



+ e l 



(13) 



t' is the renormalized hopping parameter of Eq. (|3]) . This hopping term is combined with 
a uniform gap to give in real space with sublattice co-ordinates r 



h 



3;r,r' 



M 5 , + Xr 



(14) 



for the third term in Eq. (jHj). For the diagonal term M represents a potential contribution 
of the doping atoms on the sublattice (i.e. there is a potential difference of 1M between 
the two sublattices) and t' is the contribution of the doping atoms on the next-nearest 
neighbor hopping term. After Fourier transformation we get for /i 3 



h 3 = M + x k = M + 2t' J2 cos (<V k - 



(15) 



i=l 



The term /i 3 cr 3 in H of Eq. (JSj) opens a gap in the dispersion J hi 



hi 



h\ if it is nonzero 



at the nodal points (valleys) of Jh\ + h\. Since hs is also k dependent it can vary the gap 
independently for the two valleys. It can be adjusted such that the gap parameters have 
different signs or the gap vanishes at one valley but provides a gap at the other valley. 
For instance, with = — 7r/4 and M = 3t'(l + \/3)/\/2 this term vanishes at the nodal 
points hf and hf and remains nonzero at the nodal points bj" and hf. The magnetic 
flux through a unit cell is shown in Fig. |3] and the related spectrum is presented in Fig. HI 
This model is similar to the periodic flux model introduced by Haldane [16] . However, 
the spin texture leads to some differences. In Table [1] we compare the gap-opening term 
hs of the spin-texture model with the corresponding term in Haldane's model. 
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Figure 4. Upper band of the full spectrum of the tight-binding Hamiltonian with the 
nonuniform gap defined in Eqs. (|7|), (|15[) . 



Haldane's model 



spin texture model 



h*: 



M-x'k 



X' k : 2t 2 sin0E/ B sin(k-c i ) 



2? Ef cos(k • Cj 



Table 1. Comparison of Haldane's model with our model in terms of the symmetry- 
breaking quantity /13. The vector Cj connects next- nearest neighbor sites. The 
summation runs either over both sublattices (y\- ) or only over sublattice A Q^, ). 
The definition of the variables of Haldane's model are those of Ref. [T5] . 



Below we will examine how the broken valley symmetry affects the transport 
properties. For this purpose we focus on = — 7r/4 and compare the results with gapless 
graphene (M = if — 0) and with graphene with uniform gap (M ^ 0, if — 0). 



3. Density of states 

The density of states (DOS) is the imaginary part of the single particle Green function 

1 /■ r„ „ .^\-id 2 k 



p(E) = -Im Tr 2 - / (ff-E + iO- 

7T JBZ 



n 



(16) 



BZ 



Here, the operator Tr 2 acts on the pseudospinor space and the momentum integral is 
taken over the Brillouin zone (BZ) of the honeycomb lattice, whereas Q BZ denotes its 
volume. 
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Figure 5. Density of states of the full tight-binding model is evaluated. Solid (red) 
line corresponds to gapless graphene; dashed (black) line is the DOS of the model with 



uniform gap (\h 3 



0.4i, t' = 0); dotted (blue) line represents the DOS of the 



model for nonuniform gap (£' 7^ 0) with the same m. The inset shows the same curves 
for small energies. The van Hove singularity for gapless graphene lies at the hopping 
energy E = t. The shifting of the positions of the van Hove singularity as well as the 
change of the slopes at low energies should be noticed upon a change of the parameters 
m and t' '. 



For the low-energy model of Eq. ( flOl) the DOS can be evaluated explicitly when the 
integration over the Brillouin zone is replaced by one over a circular area with radius A: 

A 



p{E) 




1 



+ 



1 



E h - »0+ E + E, 



i0- 



k 



dk 
2^2 



which gives after integration in the limit A — > 00 



p(E) = %- [&(E 

Z7T 



\m\)-Q(-E- \m\)] 



;i7) 



'ii 



This means that there are no states with nonzero energy for —\m\ < E < \m\, where m 
is the effective gap determined by the value of /13 at the gapped valley. For the case of 
gapless graphene, i.e. for M — t' — 0, Eq. (fl8|) reduces to 



p(E) = E- [Q(E) 

Z7T 



Ztx 



(19) 



The evaluation of the DOS for the model with the full spectrum and both uniform and 
nonuniform gap (i.e. with |/i 3 | = m) can be performed numerically, since the integration 
has to be carried out over the whole Brillouin zone. The DOS of graphene for different 
gap realizations is shown in Fig. |5j For small energies the DOS of gapless graphene shown 
as a (red) solid curve exhibits a linear behavior of the single Dirac cone model predicted 
by Eq. (fT9l) (note the inset in Fig. |5]). The DOS of the gapped model shown as (black) 
dashed line is zero for energies smaller than the gap m and experiences a jump from zero 
to a finite value (\m\/2ir in the linear limit) at the energy E = m. The DOS of the 
model with the nonuniform gap depicted as (blue) dotted line in Fig. [5] demonstrates an 
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intermediate behavior in comparison to the other models. On the one hand it behaves 
linearly for energies smaller than the gap m. But the slope of the linear asymptote is only 
half of the corresponding value of gapless graphene, since the number of preserved Dirac 
cones in the Brillouin zone is half the number of Dirac cones of the gapless model. On 
the other hand DOS also reveals a discontinuous behavior at E — m where the energy 
becomes sufficient to allow electrons to tunnel through the potential barrier of 2m that 
separates valence band from the conductivity band. The position of van Hove singularity 
related peaks in the DOS reveals dependence on the gap m and for m = lies at the 
next-neighbor hopping energy E — t as it can be seen from Fig. [5j 

4. Optical conductivities 

Within the linear response approach the conductivity tensor for Hamiltonian H is given 
by the Kubo formula. Then we obtain for the conductivity at inverse temperature (3 

, , i 



(Tfiv{U) 



n 



lim/ ]T (k,l\ Jfl \k,l')(k,l'\ Ju \k,l) 



y 1 f p {E kJL ,) - f p (E ktl ) d'k ^ 

Ek,i — E k)l i E k> i — E kjt i + uj — ia Vt BZ ' 

where fp(E) = 1/(1 + exp(/3(E — Ep))) the Fermi-Dirac distribution at the inverse 
temperature (3, Ep is the Fermi energy, and u is the frequency of the external field. 
\k,l) is eigenstate of H with eigenvalue E k j = (—l) l E k . The index I refers to the upper 
(/ = 0) and lower (I = 1) band, respectively. Moreover, the current operator reads 

jp = ie[H,rp] . 
For the off-diagonal matrix elements with I'^lwe obtain from the Hamiltonian in Eq. ([5]) 

(MjJM'XMUIM = Vj ^P L (E 2 k 5 ab - Kh b )H-l)He abc ^f^ 

E k v ' E k 

dh a 



(21) 



Vv,a 



dk 



The low-energy (Dirac) Hamiltonian, defined in Eq. ( TTUl) with h a = k a (a = 1, 2) and 
h% = m, gives the following expression for the current tensor: 



(kj\jn\kj'){k,l'\j u \kj) = e 






E k E k 



(22) 



In case of the low-energy approximation the integration over the Brillouin zone can be 
approximated by a circular area with radius A. The longitudinal conductivity (/x = v) 
is real due to the diagonal expression in Eq. (1221) . On the other hand, for the Hall 
conductivity, where fi ^ u, the matrix element is complex: 

(kJlj^nikJ'lj.lkJ) = e 2 ^J^, 
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<kfj,k v — (—1) ie^ u3 mE k 



(k,l\j tl \k,l')(k,l , \j u \k,l) = -e 2 



El 



10 
(23) 



From the Kubo formula of Eq. (1201) we then obtain for the real part of the conductivity 



<V(WJ 



e 



El+m 2 



[UE k ) - U-E k )} [S(u - 2E k ) + 5(u + 2E k )} kdk . (24) 



Finally, taking spin and valley degeneracy into account (factor 4) and sending A — > oo we 
obtain ES [291 



<W W J 



7re 
~2h 



(2m) 



W 



L/>(w/2) - f p (-u>/2)] [0(o; - 2m) + e(-w - 2m)] , (25) 



which can be expressed as function of three dimensionless variables as 



<W*>J 



7T6 

—f(/3m,Pfj,,Pu). 
2n 



(26) 



Hence the longitudinal optical conductivity vanishes for |a>| < 2|m| and is nonzero for 
\oo\ > 2|m|. For zero temperature the Fermi functions can be replaced by Heaviside step 
functions liing.^ fa(±u) = 0(±a;) and we obtain for the ac- conductivity the well-known 
result GS M 



Vw\u) 



ne 
~2h 



(27) 



for frequencies within the electronic bands. This value has been found to be very robust 
and barely dependent on the temperature or quality of the sample. 

One of the motivations for this work was to study the Hall conductivity in case of 
a broken valley symmetry, where we have a gapped and a gapless valley. Using the off- 
diagonal current tensor of Eq. ( 123|) . the integration over the Brillouin zone is approximated 
again by the low-energy Hamiltonian. For a uniform gap (i.e. for t! = 0) the rotational 
symmetry of the problem is restored and therefore the contribution of k^ky in the off- 
diagonal current tensor vanishes. Thus, after sending the cut-off A — > oo, we get 



/ v era f 

VM = -V37- ^ J 

M 


'ME) 


~ fp(-E) 


»( X 


1 \ dE 

+ 2E-iO+) E ' 


" yu-2E-iO+ uj 


E F \ < \m\ and T = the integral gives 


nrn 
*i» = V3 Anu lQ g 


2/ 
2/ 


n + u 
n — u 





(2? 



(29) 
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Figure 6. Real part of the longitudinal optical conductivity at half- filling for T = OK 
and m = 0.4i in units of e 2 /h. Here we show optical conductivities for: a) gapless 
graphene (red solid line); b) model with uniform gap (black dashed line); c) model 
with nonuniform gap (blue dotted line). The inset shows the same curves for small 
frequencies. Positions of the van Hove related peaks coincide with those of DOS if one 
relates frequencies to the energies from Fig. [5] by u = 2E. 



At frequency u = 2m there is a logarithmic divergence in the Hall conductivity. For large 
frequencies u this expression approaches zero, while for any nonzero gap and u — > it 
approaches the constant value that depends only on the sign of m 



m 



a 



sgn(m) e 



12 



(30) 



47r|m| h 1 h 

These are the well-known Hall plateaux of Dirac fermions ±e 2 //i |17j . This result implies 
for two valleys with gap parameters m and wl [16] 

e 2 
[sgn(m) — sgn(m')] — . (31) 



°12 



sgn(m')] — • 



5. Discussion and Conclusions 

An evaluation of the DOS (Eq. (fT6"l) ). the longitudinal conductivity and the Hall 
conductivity (Eqs. (T2U]) . (12T|) ) are performed numerically for the full tight-binding 
spectrum. The results are presented in Figs. [MZI as a function of energy (DOS) or 
as functions of frequency (conductivities) for nonzero temperatures. We compare the 
situation without spin texture it! = 0), both for the gapless case m = and the uniform 
gap m 7^ and the situation with spin texture (f ^ 0). In the latter we fix the Berry 
phase by = — 7r/4 and the uniform gap parameter by M = 3t'(l + v3)/v2 and call this 
the asymmetric valley, whereas the case t' = is the symmetric valley. 

DOS: The effect of the periodic magnetic flux, as described by h% in Eq. (I15J1 . is 
reflected by the DOS in Fig. [5j There is either a full gap (symmetric valleys) or a gapless 
DOS (asymmetric valley). The DOS of the asymmetric valley has a reduced slope in 
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Figure 7. Real part of the optical Hall conductivity at half-filling for T = OK, 30K 
and 300A' in units of e 2 /h. Numerical results calculated for the full tight-binding model 
(dots and crosses) and compared with the low-energy approximation of Eq. (|29[) (red 
solid line). 



comparison with the case where both valleys are gapless (M = t' = 0), since only one 
valley contributes at low energies. The position of the van Hove singularity depends on 
the parameters M, t' and moves to higher energies as we switch on M and t' . 

longitudinal conductivity: The behavior of the longitudinal optical conductivity in 
Fig. |6] for M = t' = reproduces the constant universal value a' = 7re 2 /2h at low 
frequencies. At higher frequencies it increases due to the van Hove singularity, just like 
the DOS. This is also the case for M^O, where the gap for t' — leads to a vanishing 
conductivity for frequencies less than the gap width A = 2m and jumps to higher values 
than the gapless conductivity at w = 2m. For the asymmetric valleys the conductivity 
behaves similarly, with the conductivity reduced by a factor 1/2 inside the gap though. 

Hall conductivity: The form of h 3 in Eq. ( TT5|) allows us to change the gap parameter 
at the two valleys separately. This means that the Hall conductivities of the two 
valleys either substract each other (for sgn(/i 3jl ) = sgn(/z 32 )) or add each other (for 
sgn(/i 3jl ) = — sgn(/z 3)2 )), when h 3 j is the gap parameter at valley j (cf. Eq. (l3Tj) ). The 
reason is that the Hall conductivities can be evaluated at each gapped valley separately, 
using the low-energy result of Eq. (130]) . Thus, for the asymmetric valleys, where one 
valley is gapped and the other is gapless, we have only a contribution of e 2 sgn(m)/2h 
from the gapped valley. This is what we see at low frequencies in Fig. [7J Moreover, 
there is a logarithmic singularity at u = 2m. It is not related to the van Hove singularity 
but appears when the frequency of the external field reaches the gap energy. Apparently, 
the optical Hall conductivity increases dramatically as the states at the edge of the gap 
start to contribute to transport [30]. The position of this singularity is quite robust and 
only determined by the gap width, in contrast to the parameter dependent position of 
the van Hove singularity. In particular, the properties of the optical Hall condutivity do 
not change over a wide range of temperatures (cf. Fig. [7]). This remarkable effect can be 
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used to measure the gap within a transport measurement. This offers an alternative to 
other methods of measuring spectral properties through transport properties [31] . 

In this work we have suggested a possible valley symmetry breaking by a periodic 
magnetic flux. The latter is generated by doping of the graphene sheet with spin 1/2 
atoms. The periodic flux opens gaps at both valleys whose values can be controlled 
independently. However, in contrast to a homogeneous magnetic field, it does not create 
Landau levels. Consequently, at most two Hall plateaux can be observed. Our calculations 
have revealed that a gap A = 2|ra| in both valleys creates the usual Hall plateaux with 
<Ti2 = e 2 /h when the gap parameter has opposite signs at the two valleys but a vanishing 
Hall conductivity if the signs of m are equal. The Hall conductivity is a^ = e 2 /2h if one 
valley is gapped and the other is gapless. Moreover, the optical Hall conductivity has a 
logarithmic singularity when the frequency of the external microwave field reaches the gap 
energy. This singularity is very robust and should be visible even at room temperature. 
Therefore, it can be used for an accurate determination of the gap by measuring the 
optical properties of graphene. 
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